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' We construct asymptotic expansions of Laplace type for the time-dependent quantum averages for 

, Bose systems with many degrees of freedom, initially populated in coherent states. These solutions 

. are localized in phase space, and they are different from the usual oscillating asymptotics for the 

' quasi-classical wave functions. These expansions are valid on any fixed time interval, and caustics 

do not contribute to the asymptotics. 
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In the present paper we develop an asymptotic theory for the partial differential equations describing the evolution 
of averages in coherent states as h — ► 0, for Bose sistems with many degrees of freedom. The crucial observation is that 
Laplace asymptotics, describing the quantities exponentially localized in the parameter of a coherent state and not the 
standard oscillating WKB asymptotics play the central role in the theory. Remarkably the phase and the amplitude in 
^vq , Laplace asymptotic expansion can be constructed using the classical Hamilton-Jacobi theory without a usual difficulty 
associated with caustics. In a sense the expansion we describe below corresponds to a "purely imaginary" WKB with 
a nonnegative phase. We use coherent states as initial wave functions of the evolution. Coherent states are the 
closest states to the classical systems because they realize the minimum uncertainty relation. They were introduced 
in R.J. Glauber Different approaches to the dynamical quasi-classical theory were reviewed in [2|, 
Our theory is based on the algebra of non-commutative operators which was discussed in detail inj 4 , 
Mathematical foundations of the theory are thoroughly developed in the works of V.P. Maslov |J, jL0|, L. Hormander 
p| and many other authors. The theory we develop below is formulated in a general form and can be applied to 
rather general Bose systems with many degrees of freedom. In particular, it can be used for computing dynamical 
quantum effects in complicated Bose-Einstein condensate systems. 
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We consider a time-independent Hamiltonian 

H(a\, ... ,a' N ,ai, .. . ,a N ) = ^^H es a)'a s , (1.1) 

1,8 

where 

e=(£ 1 ,...,£ N )€Z% , 8 = ( Sl ,...,s N ) eZf ; 
Z+ stands for the set of nonnegative integers; at* = a\ ei ■ ■ ■ a]^ N ; a s — aj 1 • • • . Here as usual 

t 1 ( t d \ 1 ( t d 

are the creation operator and the annihilation operator respectively; [at, at] — —Skifi for k,£ = 1, . . . , N. 
We assume H is formally self-adjoint, i.e., 



His = H* t 



(1.2) 
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The Heisenberg equation for an arbitrary operator valued function F(t) is 



F=~ n [H,F]. 

We introduce Poisson vectors <& Q and coherent states \a) as follows: for any a = (ai, . . . , chat) e C 



= (nh)- N/i exp^-^ (x 2 - 2V2 x ■ a + a 2 ) | 



\a) = exp 



2ft 



Then, a$ a = a$ Q , i.e., 



= 4$ Q , fc = 1 

da k 



k = 1, . . . ,N . Also 
,N . 



N 



(1.3) 

(1.4) 
(1.5) 



(1.6) 



Unlike $ a , the vectors \a) are normalized: (a \ a) — 1. Notice that $ Q is holomorphic with respect to a £ C N . 
Our goal is to study the expectation value (a\F(t)\a). To obtain the partial differential equation for this quantity we 
argue as follows. 

Let F(t) = E m , ge z* F mq (t)at m a«. Using Q, Q, ©, 



dt 



($ Q |F$ Q ) 

^($ Q |(i/F-Fi/)$ Q ) 

i £ H es F mq (t){(<£ a \a^a s ^ m a«S> a ) - (* a | a t"WV* a )} 
l - H es F mq {a* e a^a^^ a \a^ a )~a* m a s (a^^ a W^ a )} 

- H es F mq \a* e ai(h ° 



da* 



da 



J_Y( k 9_ 

"da* I \ da 



\ £ H ^ nq {a^{n^) S a^-a^n^-) a* m a q ^ exp 



da 



d_ 

da* 



In a slightly more general setting, let 



For H as in (LI) this means 



H = op Wlck H(z*,z) 



exp ■ 



(1.7) 



(1.8) 



Then from (1.5), (1.7) we get a closed equation for the evolution of (a|F|a): 

d 



dt 



(a\F\a) 



+ a ,a 



(a\F\a) 



(1.9) 
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Here, as in the theory of pseudodifferential operators, for any appropriate symbol 6(a,£) 



op[6(a,0]/ = (27r) 



-2JV 



i Rc £-(a— 7)' 



For the purpose of the present paper it will be sufficient to treat the case when Ti is a polynomial in z* , z. We will 
address the mathematical problems related to a more general TL in future publications. 
In case Tt(z*, z) is a polynomial, the equation (l.E) has the following explicit form: 



.f-V 

r! V \daJ 



H 



d \ r 



/_d \ r 
\~da~ 



H(a*,a)(f 



d_y 



(a\F\a) 



(1.10) 



Here, for r = (n, . . . , rjv) s , r! = ri! • • • r^l. 



The goal of the paper is to present an asymptotic solution to (jl.9|), ( l.lCj ) as h — > 0. We first treat the case of the 
initial condition of the following nature: 



(a\F(0)\a) =a* m a q exp 
This corresponds, i.e., to the choice 



\a - a \ : 



a E C 



A? 



F(Q) = a^ m \a )(a \a q 



(2.1) 



(2.2) 



Indeed, for a S C , 



(a|o tm |a )(ao|a 9 |a) = a* m a q \ (a , a) \ 2 
= a* m a q exp 



aga _ K| 2 _ jap 
27l 2ft 



a* m a q exp 



I a - a \ 2 



To solve asymptotically the problem (fL~9|), ( |l,10 ) with the initial condition as in (2.1) we propose the following 
procedure. Let f(a*,a,t) = (a\F(t)\a). 
Assume formally 



f(a*, a, t) = e h V^6j(a*,a, i)ft J ', 
i=o 

where the phase <S(a*, a, t) and the coefficients bj, j = 0, 1, 2, . . . are to be determined. 
The initial conditions are 

S(a*, a, 0) = — \a — ao| ; 
60 (a* , a, 0) = a*'"a 9 ; 
&,•(«*, a, 0) = for j > 1 . 



(2.3) 



(2.4) 
(2.5) 
(2.6) 



With these choices (2.1) will follow immediately. We substitute (2J3) into ( 1.10 ) and obtain the Hamilton- Jacobi 
equation for S and a system of transport equations for bj, j > 0. 
For the phase S we have from Taylor's formula: 



dS 



S = i{ H [a*, a + — - H [a* + — , a 



da 



dS 



da ' 



(2.7) 
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It is a rather remarkable fact that the Hamilton-Jacobi equation (2.7) is real (see (1.2)) and therefore the classical 
Hamilton- Jacobi theory applies. The expansion ( |2.3| ) describes Laplace asymptotics unlike WKB asymptotics valid 
for the Schrodinger equation (see, e.g., ||) describing the evolution of a wave function. 

We introduce the momenta p,p* and the effective Hamiltonian associated with (Ej) as follows: 



W(a* : a,p*,p) = -i{H(a*,a + p*)-n(a* +p : a)} 
The Hamiltonian dynamics in these variables is 

8H 



(2. 



da% 



Pk = i< 


r dn 
,dc7 k {a 


a - 


fp*)- 


dn. , 

(a - 

da k 


f-p, a)j 


.dn 


a - 


-P*) 






P*k = i< 


raw, , 


a - 


f-p*)- 


dn , , 

— a - 


hp, a) I 



(2.9) 



We now fix ao € C and solve the equation (2.7). Let 



Un .n.p .,,) = / <j [p^-( a * +p,a) - p*^-{a* ,a + p* 



n(a* +p,a) +n(a*,a +p*) 



Denote by g l the Hamiltonian dynamics corresponding to 



otk = i-^—;(o!*,a) , k = 1, 



Under the condition \n(a* , a)\ > ci\a\ M — C2, where M > 0, c\ > 0, C2 > this dynamics is defined for all a(0) € 
t G K 

a= 5 *(a(a)) 

and is a 1-parametric group of (real)-analytic diffeomorphisms ^ C^.Jl^l 

Let a belong to a small enough neighborhood of g f ao- We solve (|2.9|) with initial data 



a fe (0) , fc = l,...,JV 

« fc (Q) = -K(0) - a* ok ) , k = l,...,N 



(2.10) 



Here a(0) is chosen in a small neighborhood of ao in 

= 0, we ol 

-(a*,a)A 



Linearizing (2.9) along a trajectory with p = 0, we obtain 

d 2 n 



At 



d 2 n 



datdc&r 



datda\ 



■(a*,a)(A* r + B r ) 



(2.11) 



where A k — 5 a kl B k = Sp k , k = 1, . . . , N. 

From ( 2.10 ) the initial condition for ( [2.11 ) satisfies 

B k (0) = -At(0) 



k = 1, 



(2.12) 



We want to show that the map that assigns a{t) to each a(0) in a small neighborhood of ao where (a(t),p(t)) solves 
(2.9) with initial condition ( 2. 10] ) is a local diffeomorphism (of class C u ). Once this is established, S can be defined 
in a small neighborhood of g'ao as follows: 



(2.13) 



S(a*(t),a(t),t) = -\a(0) - a \' + / L(a*(r), a(r) , p* (r) , p(r)) dr 
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Also, p*(t) = ^pr(a*(t),a(t),t), p(t) = ^(a*(t), a(t), t). To verify the statement above we notice that fl2.1l| ) yields 



dt 



(A k B k + A* k B* k ) 



d 2 n 

datda r 



d 2 U 



■A r B k 



d 2 Tt , „ . d 2 "H . n , 

+ B r )B k + - — —A k Bl 



da* k da 



da*da k r k da r da k 



(A r 



Bt)B\ 



da k da r 

-^A*B 

datdal k r 



d 2 n 

da k da* 
d 2 H 



■A k B r 



da r da* k 



AiBt 



i~ — ~ — ~B r B k — i— — — — B*B k 



da* k da*. 
~Jt BkBl 



da r da k 



Integrating and using ( |2.12| ) we get 

A k (t)B k (t)+Al(t)BUt)+B k (t)Bt(t) 



-2|B(0)| a + |iJ(0)| s 



\B(0)f 



this implies A k (t)B k (t) + A* k (t)B*(t) < -\B(t)\ 2 . By using the second equation fl2.11p and we get C\A(0)\ > 

\B(t)\ > C~ 1 jv4(0)| with some C > uniform over any finite time interval [0,T]. Therefore, \A(t)\ > C _1 |A(0)| and 
one would never encounter a caustic moving along the trajectory (2.9) with p = (i.e., along g Oo). This is enough 
to conclude. The solution ( 2.13 ) to (2/7), (2.4) is unique in the obvious sense. 
Clearly 5((ff*ao)*, #*ao, *) = and 

j^-(( ff 'a )*, 3 W) = , ^(07*00)^*00, t) = . 
Also, the function S(a*,a,t) satisfies the inequality 

S(a*,a,t) < -Cla-g^ol 2 

in some possibly smaller neighborhood of g t ao. This follows from the above argument since A k (t)B k (t) + A k (t)B k {t) 
is the second differential of S(t) at the critical point ((ff'oo)*, <7*Oo)- 



W e now turn to the system of transport equations for the coefficients bj (a* , a, t) , j > 0. We have, substituting 
(2.3) into ( 1.10 ) and collecting terms of order h°: 



'-Ep {(!;)*<"• 



a 



/ dS\ r - 1 ? 9 
V9o" 



1 9 f/dSv-i 



< )a* e ^° + 2 da* e I (~9a 



b 



i a //9SV- 1 * 



V 9a / 9a^ 2 da( V V 9a / 



bo 



(3.1) 



In other words, 



j, v^I 9 ^* , 9S ) a ^ r ^ ^ d \ h 



N N 



=1 

9 2 



2 f-f I 

— 1 m— 1 v 



2 ,// * a s\ d 2 s 



9a* / da\da* m 



da\da* m 



U a 



95 \ 9 2 5 
9a ' / dagda r 



bo 



= Lob, 



(3.2) 



In the right side of (p/fl) It, I = 1, . . . , N is the vector with components 5^, s = 1, . . . ,N. 
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For an arbitrary j > 



bj = Labj + ^ L kbj-k ■ 



(3.3) 



fc=i 



The differential operators in the right hand side of ( |3.3|) are as follows: 

i d k+1 ^ 1 f/S 



-1 r \ 



{k + l)\dh k + x r! I \da 

re 



— H(a,a){h- + - 



(3.4) 



ft=0 



To solve the system (3.1)— (3.4) with initial conditions ( p.5[ ), (2.6) we use (2.9) to obtain with fixed a*(0), a(0) in a 
small neighborhood of a 



-&o(«*(*), «*(*), t) + co(a*(t), a(t), t)b = 
at 



(3.5) 



where cq is the coefficient in the second term of ( |3.2| ). Having solved (3.5) with initial condition ( |2.5| ), we obtain a 
first order linear inhomogeneous ODE for b%, etc.. According to ( |3.3| ) to solve the equation for 6j we are recursively 
faced with the problem of solving the same ODE with a different known right hand side. The solution to this problem 
therefore exists and is unique for (a*, a) being sufficiently close to ((g t do)* , g f oto) ■ 
To solve the problem (1.9), fll-lCP with the initial condition 

(a\F(0)\a) =a* m a q 

we use the completeness relation 

^ |ao) (a | d 2N a = id 



(TTft) 



which leads to the expansion 



(a\F(t)\a) = 



(nh) 



N 



d 2N aoe i°" t} bj(a*,a, t)K 

3=0 



Here in the right hand side the phase S and the coefficients bj , j > Q depend on a parameter ao G 



4. Example (see [12], [p| for the exact solution) 



Let N = 1, Tl(a* , a) = uja*a + fia of. The equation (2/7) takes the form 



dS OS 



The first transport equation is 



, • / * 9 d \ , _. / ,JS 9 2 dS d . 



d 2 S 



,d 2 s 



b = L b 



da* 2 da 2 
A simple computation yields: 

S(a*,a, t) = -\a(0) ~ a a \ 2 - i\xt (a*(0) 2 o§ - a{0) 2 a* 2 ) 

+ |«(0)| 2 ^1 — e - 2 Vt(a*(0)a ~"(0)aS)^ 

b (a*,a,t) = a(0)* m a(0)«- 



exp 



i(a(0)*ao — a(0)o;g 



2|a(0)||oo| 



arctan(2/if|ao| |a(0)|) > , 



(4.1) 



(4.2) 

(4.3) 
(4.4) 



where 

a = a(0) exp{iwi + 2nta* o a(0)} . (4.5) 



Notice that for any t the map (4.5) a(0) — > a is invertible in a sufficiently small neighborhood of ao 



Conclusion 

In conclusion, we developed a theory for Bose systems with many degrees of freedom which describes not an 
oscillating asymptotic behavior of a quasi-classical wave function, but a Laplace asymptotics for quantum average 
quantities, which are localized in phase space. This approach has many advantages in comparison with the standard 
one, as it does not require (i) an intermideate step of computing a complicated behaviour of a quasi-classical wave 
function which usually presents difficulties associated with caustics, and (ii) an additinal integration for computing 
the average values. Our approach can be used for solving many dynamical problems related to Bose systems with 
many degrees of freedom. 
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